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1. INTRODUCTION 


The Agriculture and Resources Inventory Surveys Through Aerospace Remote 
Sensing (AgRISTARS) program Is a major effort to assemble and apply the remote 
sensing technology to agricultural analyses and Earth resources Inventories. 
AgRISTARS does not rely on remote sensing alone; It also Integrates available 
ground observations and meteorological data to provide an enhanced perspective 
of crop conditions ai.d Increased capability for monitor ng temporal and spatial 
changes, "he meteorological conditions and their Impact on the environment are 
vital to the detection of the agricultural conditions. Weather data, there- 
fore, become a major component of the ground-truth Information required for the 
effective use of remote sensing technology. 

The Foreign Commodity Production Forecasting (FCPF) project within AgRISTARS Is 
particularly concerned with the application of remote sensing technology to 
agriculture In foreign areas. The FCPF effort Is directed toward Inventorying 
and monitoring the production of crops that are significant to the world' ? food 
supply. In few applications does the role of weather have such a tremendous 
Impact as In agriculture. The project, in response to this fact, has required 
weather analysis and Interpretation activity. This activity has specific 
objectives In support of FCPF experiments and technology. 

The primary objective of this task is to provide relevant information about 
meteorological conditions In all project experiments or similarity regions. 

This includes an assessment of the impact of weather on crop development, 
yields, and spectral appearance. Task Input relies upon the production of 
meteorological summaries and ui splays of various meteorological variables and 
Indices. These products aid In the Interpretation and assessment function. 

A supporting task objective Is to access, develop, and manipulate all available 
meteorological data. A computer capable of handling large volumes of data is 
required. The attendant software must be compatible with the mode and format 
of varied data sources. It also must be flexible in order to respond to 
different requirements for Information. 
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A third objective Is to provide date and Interpretive Information of quality 
and detail. This requires procedures to correct Inconsistencies In the data 
and to fill any data hiatus detected. This Is accomplished with software based 
upon operational experience w1*h the characteristics of meteorological observa- 
tions and knowledge of the vagaries of the spatial distribution of the data. 

The final objective is to provide a set of output products which meets user 
requirements; employs pertinent parameters; and displays all features In a 
standard, unambiguous manner. These products will be the basic tools of the 
overall agricultural Interpretation and assessment activity. 

This report provides a description of the current line of products which was 
devised and delivered for the U.S. /Canada Wheat and Barley Exploratory 
Experiment. The functions and/or applications of these products are also 
given. The data requirements to produce these Items are exhaustive and are 
described in detail. Procedures for assembling and processing these data are 
also presented. The software to accomplish this task, which Is still being 
developed and streamlined. Is documented In Its current state. 
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2. PRODUCTS AND FUNCTIONS 


The primary purpose of meteorological data Is to provide point-specific Infor- 
mation, whenever possible# Weather analysis and Interpretation procedures 
developed during the Large Area Crop Inventory Experiment (LACIE) were regional 
and relied on interpretation of network weather observations# The approach for 
the AgRISTARS FCPF project Is to provide segment-level detail and analysis of 
weather data In a time series. This method eliminates large volumes of 
extraneous regional Information requiring subjective Interpolation and replaces 
It with segment-specific data obtained objectively. Temporal variations are 
reduced to a weekly time series. 

To display the desired data for a particular site over the growing season, a 
three-part cllmagraph Is employea. A Fortran program Is used to generate this 
product on a line printer. One cllmagraph Is produced for each segment 
location In the experiment region. 

Part 1 of the cllmagraph consists of a composite plot of (1) weekly observed 
and mean temperatures and (2) weekly and normal precipitation totals (see 
fig. 2-1). The first header line specifies the data acquisition year, the 
segment number, the state and crop reporting district (CRD) In wh1c* the 
segment Is located, and the latitude and longitude of the site. Ths second 
header line lists the available Landsat acquisitions. The horizontal time 
scale Is the Julian date for the ending day of the week. The temperature scale 
Is degrees Fahrenheit at 2° intervals, and the precipitation scale Is In Inches 
to the nearest two-tenths of an Inch. Climatological normals are plotted as 
asterisks (*); observed values, as pluses (+) for above normal and minuses (-) 
for below normal values. These four parameters are not observations or records 
at the segment location but are Interpolations of values from the available 
first -order network for the region. The values at the stations are summarized 
by week prior to interpolations. For this particular plot, 26 Interpolations 
were needed on 4 variables to provide the segment-level estimates. 
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Figure 2-1.- Climagraph product, part 1 


The cllmagraph display of data penults Interpretation of potential Interactions 
of temperature and precipitation with spectral appearance. For example, the 
week prior to acquisition date 175 was cooler and wetter than normal. This 
provided some relief from any ongoing stress, as well as conditions for a 
reduction In the rate of cr:p development. 

The plot also aids In detection of long-term effects and sequences of signifi- 
cant weather features. For example, for 5 weeks prior to day 238, temperatures 
were consistently below normal. This phenomenon probably delayed ripening of 
spring wheat. The heavy rainfall during the week ending on d:y 210, an extreme 
departure from normal, would probably have been detected as standing water In 
some fields If an acquisition had been available for that date. 

Part 2 of the cllmagraph depicts the crop moisture Index (CMI) and the adjust- 
able crop calendar (ACC) model outputs by week (see fig. 2-2). The header 
Information, available acquisitions, and weekly time scale are the same as In 
part 1. The scale for the CMI Is dimensionless with a normal reference line at 
zero. The ACC scale. In this case. Is for the Robertson blometeoro logical time 
scale for spring wheat (~ef. 1). The stages In this scale are: planting » 1.0; 
emergence * 2.0; jointing > 3.0; heading * 4.0; soft dough • 5.0; and ripe • 6.0. 
This scale Is also dimensionless and Is used to depict the cumulative develop- 
ment of the crop as estimated In a meteorologically driven model. Again, the 
26 values of each variable are interpolated estimates derived from running the 
models with the available first-order station data. 

These two plots characterize the response to the meteorological conditions at 
the segment level and facilitate the assessment of meteorological Impacts on 
spectral appearance. The CMI Integrates the effects of both temperature and 
precipitation Into a budgeted Index related to soil moisture and, hence, 
stress. For Interpretation of the Index values, consult table 2-1. 

The CMI plot when compared to the temperature and precipitation plots Indicates 
to some extent the response of the CMI to the other two variables. For 
example, the week ending on dqy 168 was warm and dry enough to begin slight 
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[As reported in the Weekly Weather and Crop Bulletins prepared 
by the U.S. Department of Commerce, Department of Agriculture] 
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Figure 2-2.- Climagraph product, part 2. 


stress* but conditions recovered In the subsequent cool* wet week ending on 
day 175. In general, the Index did not Indicate stress, despite below normal 
precipitation, until the higher temperatures of mid to late summer began to 
take effect. The sequences in this time series are also harbingers of spectral 
appearance. Sustained and unrelenting stress conditions are Indicated after 
the week ending on day 224. 

The final plot of ACC values by week Indicates at what particular stage of crop 
development favorable or unfavorable meteorological episodes or sequences 
occurred. It Indicates the best current estimate of a crop's stage of develop- 
ment as derived from the current year's weather. The plot, which defines the 
growing season for the crop. Is used to establish the Intervals (biowindows) 
between spectrally significant crop stages. Acquisitions during these 
Intervals are selected for classification procedures. 

Part 3 of the cllmagraph product provides information which cannot be conven- 
iently plotted and a special synopsis of special features of the ACC plot (see 
fig. 2-3). The header and available acquisition date Information is followed 
by a table. entitled "Date and Region of Reported Significant Ecological 
Events." An effort has been made to collect reports of events or situations, 
not necessarily weather related, which may effect spectral appearance in the 
experiment region. This type of information is dated and assigned to the 
appropriate level of geographic detail. For example, localized hail damage In 
North Dakota, CRD 3, was reported for the week ending on day 252. Although 
segment 1387 Is located In this CRD, hall may not have affected the segment 
itself. Likewise, the report from the week ending on day 147 indicates that 
wild oats was becoming a statewide problem, although not necessarily in 
segment 1387. Such reports are provided to the analyst as potential clues for 
explaining anything unusual observed In the segment Imagery or spectral 
sequences. 

The next feature c* part 3 Is entitled "Adjustable Crop Calendar Dates for 
Spring Wheit by Acqu.sltion Date and Stage." These are data points which may 
be extracted from the ACC plot but which are displayed here for convenience. 
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Figure 2-3.- Climagraph product, part 3. 


The model stage estimates for each available acquisition date are given, 
followed by the estimated dates of the six cardinal stages of the crop model. 

The final Item on the cllmagraph Is the specification of potential biowindows 
for acquisition selection. In this case, two are defined in terms of the ACC 
cardinal stage dates. Window 1 Is the 24-day period from 5 days before to 
18 days after spring wheat planting. Window 2 Is the 21-day period from 
10 days before to 10 days after heading. Additional blowlndows were defined 
for barley, the other crop of Interest, but ACC models are not yet available. 

These windows are determined by comparing the historical average of the bio- 
window definition date (e.g., barley turning date minus spring wheat heading 
date) to the available ACC date. For example, if historically for the region, 
barley turning (key date of the biowindow) occurs 10 days prior to spring wheat 
heading, the best estimate In the current year for this date is 10 days before 
the ACC-estimated date of spring wheat heading. 

This explains the current version of the climagraph product. Modifications are 
anticipated in response to user requirements, to analyst evaluation of the use- 
fulness of the displayed variable, and to the introduction of new variables and 
indicators of spectral appearance. The central theme remains the careful 
reduction of data to the segment level for a specific time Interval. 

Two additional products are generated from the weather analysis and Interpreta- 
tion activity. Both contain special Interpretations of the CMI, the primary 
variable for study of crop spectral appearance. 

The regional CMI Interpretation (fig. 2-4) provides an expanded view of the 
extent of stress conditions (as Indicated by the Index) over time and space. 

The product Is a printed narrative for each subregion (CRD) In the experiment 
region with the segment numbers Included. 

The second CMI product (table 2-2) provides a subjective assessment of the 
Impact of CMI variations on crop development and spectral appearance. The 
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TABLE 2-2.- CMI INTERPRETATION FOR CROP DEVELOPMENT AND SPECTRAL APPEARANCE ASSESSMENT 
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AND F A VUx ABLE MOISTURE CONDITIONS AS INTERPRETED FROM The Chi. 
THE EASTERN CRO'b E xpEh 1 EhCED LITTLE OH NU a T REsS CONDITIONS 
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WERE QUITE DhY BY eAWLY JULY. ThEheaFTeW THINGS IMPROVED 
and CONDITIONS BECAME VERY MOIST BY MIO-SEPTEmbEW. 

1.2.7 SOUTHEAST (SEGMENTS: 1 3*9 . 1 sbB . 1 6b 1 . 1 t>6* . 1 92 a . 197 A ) 

The Cm I wAS FmIhly STablE wITm CONDITIONS AoEOUATe UW 
FAVOHABLE ThhuOGH M I 0-SePT£MBEW . 

1.3 MINNESOTA 

1.3.1 NOHTH«EST (SEGMENTS: 1S14. 151B. 1825. 19b7) 

CONDITIONS wEhE adEouaTe Oh FavOhAbLe OVER THIS REGION 
ThhOUGh JULY pITH STEAJY DETEHIOKAT Ion IN AUGUST AND 
SEPTEMBER. 

1.3.2 wEbT CENThal (SEGMENTS: 18JS. 18*21 

FAVOHABLE MOISTUHE CONDITIONS PERSISTED THROUGHOUT THE 
GHOhING SEASUn EXCEPT FOP EXCESSIVE mOISTukE IN MlD 
TO LATE JUNE. 

1.3.3 CENThal (SEv»menTS: 152*. 18*3) 

MUISTUHE PAS auEUUATE To EXCESSIVE THROUGHOUT Th£ SEASON 
PITH EXTREMELY w£T CONDITIONS PHIOH TO ThE DELAYED 
PLANTING. 

1.3.4 SOuTh. EST (SEGMENTS: 13»0) 

muISTUhE «AS AUeQUATE To EXCESSIVE ThhOUGhOUT ThE 
se ason w I Th EXT Ht.Mfc.LY -»tT CONUT IONS COMMENCING In AUGUST 
WITH LITTLE ImpwOVEMENT ThROUohOUT Th£ HARVEST PEHIOU AND 
INTO LATE SEPTEMBER. 

1.4 SOUTH DAisOTA 

1.4.1 NORTH central (SEGMENT: 1599) 

MOISTURE CONDITIONS p£nt ADEOUATE THROUGHOUT ThE 
SEASON EXCEPT FUR SOME STRESS DURING The FIRST HALF 
Of JUNE . 

1.4.2 NORTHEAST (SEumEnT: I960! 

THE CMI INOItAlED FAVURmBLE MOISTURE CONDITIONS ThhOUGhOuT 
The SEASON BUT MOISTURE -AS ExCESIVE DURING ThE LATTtR 
HALF OF JUNE. 

1.4.3 CENTRAL (SEGMENTS: I6a9.176s) 

CONDITIONS wEPt ADEQUATE ThROUGh MID-SEPTEMBER EXCEPT 
FUh SOME SLIohI STkESS IN EAhLY JUNt. 

1.4.4 SOuTh CENThal (Segment: i « s 7 e> ) 

The Cmi mas ADEOUATE an.) STablE ThhOUGh LATE SEPTEMBER. 

1.4.5 SOUTHEAST (StuMENTS: 1M*) 

MOISTURE PAS huRE Than mUEQUATE ThROUGh ALMOST ALL Tht 
SEASON ■ I T m ExE CSSIVE muISTURE INOICaTED In mID-MAY. 

1.4 MONTANA - CENTRAL ( St GM£N T : H*H) 

MO 1ST URt MAS AOEoUA T t Through mat but DtCLINED ANO REMAINED 
MARGINAL .TO POOR ThPOUGhOUT Th£ hEM a I NDt R uF Th£ SEASON 
WEACHINO NtARLY OWUUGH T CONQIIIONb oY LAft AUuUST. 


Figure 2-4.- Regional CMI interpretation. 
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assessment Is made at the specific CRD level In terms of expected deviations In 
stage and color at the six normal stage dates. This product represents a 
detailed Interpretation by an experienced analyst of appearance and the reasons 
for an anticipated spectral sequence of spring wheat In the CRD. 

All of these products were prepared and delivered In support of the U.S. /Canada 
Wheat and Barley Exploratory Experiment. They are still under evaluation and 
subject to modification. For now, these products represent the current state 
of weather analysis and Interpretation techniques In support of FCPF experi- 
ments. They can be produced for almost any region for which data are avail- 
able. These data requirements are not entirely inflexible; however, timely 
acquisition of meteorological data Is the most crucial component and must be 
supplied without compromise. 
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3. DATA REQUIREMENTS 


Data requirements fall Into two categories: (1) ancillary data and (2) meteor- 

ological data. The ancillary data are used to define a suitable subset of the 
meteorological data spatially and temporally and to control the structure of 
the Information displayed In the products. The meteorological data are 
refined, summarized, used to run models, and displayed for Interpretation. 

3.1 ANCILLARY DATA 

Ancillary data must be obtained for each point (segment) of Interest to the 
user. The types, functions, and sources of these data are listed In table 3-1. 
These data sets, which are In card Image format, do not require a large amount 
of storage space. Usually, all the data of a particular type can be combined 
for all segments In a geographic region. For example, items 1 through 5 In 
table 3-1 can be placed on a single record, and the file of these records can 
comprise all segments in a given experiment region. 

3.2 METEOROLOGICAL OATA 

Two classes of meteorological data are required by these procedures: (1) his- 

torical records of monthly data and (2) daily records from the weather stations 
for the year of Interest. 

The historical records must be a minimum of 10 years in length, preferably 
30 years. They must be at th« lowest level of representation available 
(station, whenever possible) and Include monthly averages of mean temperature 
and monthly totals of precipitation. These data are essential for proper 
Implementation of the CMI algorithm. The normals of temperature and precipita- 
tion, Items 7 and 8 of the ancillary data (table 3-1), may also be determined 
from these records. 

The dally records must include the maximum and minimum temperatures and the 
total precipitation. The records must be complete for each Item for each day 
of the period of interest and Include as many stations as possible within the 
region bounds. 


3-1 


TABLE 3-1.- ANCILLARY DATA REQUIREMENTS 
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a Header refers to identification information in computer-generated products. 




Because of the enormous quantity of historical and dally data, tape and/or disk 
storage Is required. For completeness of the records and the dally observation 
reports, further processing of raw data Is required. Ideally, the historical 
records will correspond to the dally records station by station. 

A summary of the meteorological data requirements for these procedures appears 
In table 3*2. The International scope of anticipated FCPF experiments requires 
worldwide coverage for meteorological data. The minimum station density 
requirements for these procedures have not been determined; however. It does 
not seem likely that such requirements will ever be exceeded In foreign areas. 
Therefore, observations from every available station are needed. 


TABLE 3-2.- METEOROLOGICAL DATA REQUIREMENTS 


Type 

Function and/or 
application 

Data source 

Mean monthly 
temperature 
(120 months) 

CMI operations 

Historical records and 
published values 

Total monthly 
precipitation 
(120 months) 

CMI operations 

Historical records and 
published values 

Dally maximum 
temperature 

ACC models, computing 
weekly means, and CMI 
operations 

Dally records from NOAA® 
tapes and published values 

Daily minimum 
temperature 

ACC models, computing 
weekly means, and CMI 
operations 

Dally records from NOAA 
tapes and published values 

Daily precipitation 
totals 

Computing weekly totals and 
CMI operations 

Dally records from NOAA 
tapes and published values 

Latitude and 
longitude of first- 
order stations 

Objective analysis 

* 

Station records 


^National Oceanic and Atmospheric Administration. 
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4. PROCEDURES 


An Idealized configuration for these procedures appears In figure 4-1. A user 
provides basic Inputs about the desired region through an Interactive terminal 
of the Laboratory for Applications of Remote Sensing (LARS) computer at Purdue 
University. This computer then accesses Installed data bases and software to 
obtain a machine-generated tool of analysis (e.g., cllmagraph) based on 
weather conditions for each segment. The computer system provides permanent 
disk storage, data tape handling facilities, and working temporary storage. 

The software provides extraction modes for all available data, quality checks 
of the data, summaries, and algorithm outputs of meteorologically driven 
models* 

4.1 ASSEMBLY OF INPUTS 

Prior to the creation of products at the terminal, numerous small data sets 
and ancillary Information relating to the experiment region and Its segments 
must be assembled. The first requirement Is a geographic definition of the 
experiment region. If the region has distinct geopolitical boundaries and 
does not exceed the size of a square 10° latitude by 10° longitude, the 
definition becomes the smallest quadrangle, bounded by whole-degree lines of 
latitude and longitude, which encompasses the region. 

This constraint on size has two supporting arguments. First, an area that Is 
much larger would probably not be cllmatologlcally homogeneous. The second 
reason for restricting the geographic size Is the result of adaptation of 
software to objective analysis techniques. The number of grid cells to be 
analyzed for a particular grid size Increases with area. This, In turn, 
Increases computer core requirements and computation time. If the political 
region exceeds the 10° latitude by 10° longitude size, two or more quadrangles 
must be specified. Precise Information on the climatic regime, topography, 
and spatial distribution of the segments Is required to define the size and 
positions of these subregions. The distribution of the segments within the 
region becomes a key feature when geopolitical boundaries are not definitive. 
The size constraints also apply to the special quadrangle to encompass this 
distribution. 
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LARS SYSTEM 



Figure 4-1.- Idealized configuration for weather analysis and interpretation procedures. 



The distribution of the segments within a region Is best determined from a 
complete list of the segments, along with their latitudes and longitudes, to 
be used In the experiment region. Ultimately, the coordinates will be used In 
an objective analysis to obtain estimates of various meteorological parameters 
at the segment level. 

With tiie definition of each experiment region quadrangle, the specific 
regional data must be assembled In disk storage. These data Include dally 
meteorological data for all available stations, monthly climatological data, 
estimated planting dates for each crop of Interest at each regional segment, 
estimated regional soil moisture capacities, significant ecological events, 
and open and close dates for data collection. 

The period of data collection Is largely determined from the historical 
regional crop calendars. This meteorological data collection period should 
overlap the Imagery collection period during the experiment. The dally 
meteorological data and the data for significant ecological events are 
acquired as specified by the open and close dates and may be assembled as they 
become available. 

Dally meteorological data are extracted as specified by the latitudes and lon- 
gitudes of the quadrangle and the collection period. These raw data require 
some quality control because of Inconsistencies, such as maximum temperature 
reported as less than or equal to mlnlimin temperature. Such reports are dis- 
carded and regarded as missing. All missing values are replaced with esti- 
mates, when feasible. If the observation record for a particular station is 
less than 25 percent complete, it Is deleted from the data set because values 
for too many of Its observations would have to be estimated. 

The completion of the remaining observation matrix (1 stations by j days by 
3 variables) Is attained through an objective analysis technique. A value for 
a missing observation of a particular variable on a given day Is estimated 
from the scalar field of that variable. The software can be used to determine 
objectively a unique solution of the scalar field. This solution, which Is 
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obtained from grid specifications and available observations within the 9 r 1 d t 
can then be used to estimate a value for a;jy given point In the field. This 
procedure Is repeated until a complete data set of maximum and minimum 
temperatures and total precipitation values Is obtained and stored on disk. 

Tne historical records of monthly climatological data acquired within the 
experiment region quadrangle must be assembled. If the record 1$ for a divi- 
sion rather than an Individual station, approximate geographic coordinates of 
the division must be obtained because the data must be treated like those of a 
station for computational purposes. 

Estimates of soil moisture capacities for each station within the experiment 
region must be made. This can be accomplished by consulting either the pub- 
lished data for the region or a soils expert capable of utilizing general 
soils classification maps for the region. The estimates are to the nearest 
Inch and frequently are uniform over entire regions. 

A crucial element of data to be assembled is estimates of regional planting 
dates for each crop of Interest. This may be accomplished using any of three 
resources: (1) historical normals; (2) starter models, such as that developed 

by Feyerherm for spring wheat; and (3) ground-truth data from the segments. 
From whichever means Is appropriate, dates must be estimated for the first- 
order stations. 

The final data type to be assembled prior to processing is tfte significant 
ecological events for the region. These are terse, verbal comments extracted 
from published reports about crop and weather conditions In the region. Each 
event must be Identified by the approximate Julian date of the reported occur- 
rence. The event must further be Identified by region and subregion codes; 
for example. North Dakota (ND) and northeast CRD (03). 

All the preceding data (whether compiled as cards, tapes, or disk files) must 
be assembled collectively as a set of card Images In permanent storage where 
they are accessible, during program execution. The format specification of 
these files must be consistent with ci . rent software. 
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4.2 DATA PROCESSING 


All variables that must be reduced to the segment level through objective 
analysis must be processed and available at the first-order stations prior to 
Interpolation. 1)1$ assures that all models and statistical procedures are 
applied to data with only observational errors. If the raw data for observed 
temperature and precipitation were reduced and then processed at the segment 
level, the inherent error of Interpolation would be compounded by the error of 
observation prior to running the models or producing statistics. The Impact 
of this compound error is not known, which complicates any assessment of error 
statistics of segment-level variables. 

The first stage of the processing is to obtain the daily stage of development 
for each crop of interest by running the ACC model (s) from the estimated 
planting date and the daily meteorological data for each station. One further 
step is required prior to interpolation: extract only those values for the 
week-ending Julian dates for the experiment time period. 

The next variables for processing are obtalnen by contrasting the daily 
temperatures and precipitation values to the mean temperature and the total 
precipitation for each week-ending period. These weekly values not only will 
be interpolated but also will be used to run the CMI algorithm. 

Next, the historical records are processed. Two products are obtained: 

(1) coefficients for the CMI algorithm and (2) normals of temperature and pre- 
cipitation. Both are outputs of the Palmer hydrological accounting system 
(ref. 2) operating on the entire record. 

The normals of temperature and precipitation are then utilized to obtain simu- 
lated weekly values corresponding to those determined for the current year. 

The simulation is accomplished by harmonic analysis. 

The final variable to be obtained at the first-order stations is the weekly 
CMI values. These values are computed by applying the CMI algorithm to the 
derived coefficients for the station and the weekly values of temperature and 
precipitation. 
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The following weekly variables by station are now ready for objective analysis 
to obtain the segment-level estimates displayed In the cllmagraph products: 

(1) mean temperatures, (2) normal mean temperature, (3) total precipitation, 
(4) normal total precipitation, (5) final ACC stage estimates, and (6) CMI 
values. The objective analysis of these six variables for each week utilizes 
the latitudes and longitudes of both stations and segments and the grid for 
experiment region quadrilateral. 

The final data processing at the segment level applies ancillary data and 
meteorological data to produce segment-specific cllmagraph products. These 
products are printed for delivery In support of FCPF experiments and then are 
stored on magnetic tape, along with all their contributing data files. 

4.3 INTERPRETATION OF THE CMI 

Both CMI Interpretation products are subjective in nature. The regional 
narrative Is based upon examination of CMI maps and plots to obtain a coarse 
assessment over time and space. The interpretation for stress and the impact 
on crop development and signature response are based upon an analyst's 
experience with the spectral sequence of crops. This is complemented by the 
analyst tracking the CMI values throughout the season to detect possible 
deviations from normal development In the region and to anticipate the 
spectral appearance. 

4.4 MODELS AND ALGORITHMS 

The important models and algorithms used In these procedures are listed in 
table 4-1. Fortran Implementation of these items has been completed. They 
exist as flexible subroutines accessible by an evolving line of Conversational 
Monitor System (CMS) executive processors and other Fortran programs which 
access, manipulate, and display the meteorological data in the cllmagraphs. 

The software of these models and algorithms represents the stable core of cur- 
rent Fortran programs used in these procedures. Documented program listings 
are given In appendix A of this document. The structure for Interactions of 
this software appears as a flowchart In appendix B. 
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TABLE 4-1.- MODELS AND ALGORITHMS APPLIED 


Name 

Function and/or 
application 

Data source 

Wagner variational analysis 
technique 

Objective analysis of 
meteorological variables 

Ref. 3 

Crop moisture Index 

Prime stress Indicator 

Ref. 4 

Palmer hydrological accounting 

Local coefficients for CMI 

Ref. 2 

Thornthwalte model for potential 
evapotranspiration (PE) 

PE values for Palmer and 
CMI algorithms 

Ref. 5 

Robertson biometeorological time 
scale model for spring wheat 

ACC stage estimates 

Ref. 1 


The Wagner variational analysis technique (ref. 3) is applied for objective 
analyses of meteorological scalar fields In regions of sparse data. This 
method utilizes a low-pass filter to provide a consistent and computationally 
rapid means of estimating values at grid points in the analyzed field. The 
errors associated with estimates of the variables used in the climagraphs are 
being investigated. 

The CMI is based upon a relatively simple two-layer soil moisture budget 
(ref. 4). On a weekly basis, the value of the index is computed using the 
total precipitation and an "appropriate" value for evapotranspiration (ET). 
This ET value is an adjusted value of PE computed from mean temperature. 
Coefficients for PE and the adjustment factor are based upon the Palmer hydro- 
logical accounting approach (ref. 2). This technique is applied to the entire 
monthly climatological record and assumes that an estimated soil moisture 
capacity for the region can be used in the soil moisture budget. Surpluses 
and deficits of moisture for the Index are therefore linked to normal 
conditions for the region. The method for computing PE is based upon a model 
proposed by Thornthwalte in his efforts to simplify climatic classification 
(ref. 5). 
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The Robertson blometeorological time scale model for spring wheat (ref. 1) 
quantifies the progress of the crop toward maturity as functions of dally 
maximum and minimum temperatures and daylength. The computations result 1r. a 
dally Increment of development during the crop season. These Increments are 
accumulated from an Initial value of 1.0 until a value of 6.0 (ripe) is 
reached. The six stages described earlier require five sets of coefficients 
(one for each interval). As the accumulated value reaches a new stage, the 
coefficients applied to the dally variables change, reflecting the changes In 
the response of the plant to Its environment as it matures. The model was 
tested and used successfully and extensively during LACIE for both weather 
interpretation and advanced yield modeling. 

4.5 SCHEDULING 

The procedures described above are applied to two main types of experiments: 
(1) a crop-year experiment in which the region and time period are specified 
prior to planting and the data collection and processing progress with the 
growing season and (2) a historical -year experiment in which the region and 
time period are specified after harvest and all data have been collected. A 
generalized schedule for providing support to these two experiment types 
appears in figure 4-2. 

The crop-year experiment schedule runs for about 12 months, approximately 9 of 
which are devoted to data collection during the growing season. Final 
processing begins after 11 months when the final meteorological data become 
available. At this point, coordinates for the segments in the region are 
essential to the objective analysis of meteorological variables and the 
production of the final climagraph product. 

The schedule for a historical -year experiment runs about 4 months after the 
selection of the region and crop year. Approximately 6 weeks is anticipated 
to procure the meteorological data required for the experiment. An additional 
9 weeks is required to extract the station data, run the models, receive the 
segment coordinates, perform the objective analysis, and produce the products. 
This schedule would probably have to be extended to 6 months if multiple years 
for a particular region are included in the experiment design. 
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Perform objective analysis 
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Figure 4-2.- Generalized weather analysis and interpretation 

procedure schedule. 



5. CONCLUDING REMARKS 


The procedures and products developed for the weather analysis and Interpreta- 
tion activity supporting the U.S. /Canada Wheat and Barley Exploratory Experi- 
ment are described In this report. Innovative approaches have been applied In 
order to meet the task objectives, and an Improved line of products has been 
created. These products are tools In this effort and are subject to evalua- 
tion and revision as user requirements are better defined and understood. 

The major emphasis Is to develop segment-specific Information about weather 
conditions and an objective, machine-oriented representation of this 
information. 

The data requirements of these procedures have not Increased significantly 
since LACIE; however, the requirements of accessibility and geographic 
coverage have changed. Data tape libraries and disk data bases are needed to 
assure timely processing and to support the massive data reduction described. 

The procedures will continue to change as new requirements and techniques 
evolve and as new products are designed. As data bases become comprehensive 
and stabilize In format and reliability, the software for their manipulation 
and reduction will become more uniform and streamlined. The basic software 
for models and algorithms will expand as new ones are developed by supporting 
research. Future changes probably will require revision of these procedures. 
This document represents the current state of development for weather analysis 
and interpretation techniques for the FCPF project of the AgRISTARS program. 
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APPENDIX A 


PROGRAM LISTINGS OF MODEL AND ALGORITHM SOFTWARE 


APPENDIX A 


PROGRAM LISTINGS OF MODEL AND ALGORITHM SOFTWARE 

Listings for the WAGNER, CMISUBS, and ROBBMTS computer programs are given In 
this appendix* 


A-l 


FTLF! waGNEP FOPTRAN a MUPOUF / LAOS 3031 


1MFRFO IOATa ( 1 00 »M 

0 TENSION »L»T (*"") . *LONS(GO 0 ) «X(in 0 .fl).YL»T( 500 ) . YLONG(SOO) 
niMJ:N«I^N iMSAtAftlfiUSHtSAIilfACiPliHAISl 
PFftL MAXI A T « M 1 v l AT.MA *l.M 3 ,MlnL (u 

OF AO (S. 10 1 .FM'ISI 3 ) *A»L AT .MTNL XT • M A*LN r nMj NLNG • l S I ZE . JS I ZE » L AH 
JO) F 0 P-AT(F?.fi.F 3 # 0 .?F*. 0 , 2 I 3 . IX.S 44 ) 

IF (MAll AT. lT. -40.0) ST<)K» 

GO TO ^14 
ST0» 

CONT INUF 

wo i TF ( *i . oi) l ) *MILAT«MINI AT.maxlnG.mINLNG. IST7E. JMZF 
FORMAT ( 1 M 1 . 41 . 2 ?H*M* IMIJM LATITU-F .FS. 2 /SX. 

1»?H»*INI m iim LATITill.c .F*.2 »/*x.??MMAXImum LONGITUDE • 

?F 6 . 2 /*» . 2 ?* M TNT , 'Mi , >‘ LON 3 IT IOF ,FS.?/*>X. 

3 ? 2 m^U'<ofo n-s O&n °OINTS. | 6 /«*.? 2 mnummER E-w r.oio POINTS* 16 ///. 

*• '.ISA4/////1 

ma»LAT 3 MAXIMUM I eTt TOO? ON 4 A» 

M f NL A T 3 MINIMUM L A T I P !HF ON MAP 
Mf XL NO 3 MaxTHIM LONG* TljOE on * 4 P 
M I NLNGsm jxj t mijm lOMr.TTOOF ON 
1ST 7rsN'|MMFP DC OMin in n-S OIBFCTI'XN 
j«T 7 F * NUmoFR OF f,PIO Tn E-. 0 1 RFC T I ON 


*13 
ST A 

*«1 


T S I 7 F • M A XL 4 T .mi.'Lat mijsT mc at JUSTFO SO 
fio SO 


ISC»LF 

JSCALF 


II 


an INTEGER 

AN INTEGER 


JOI7F,maxln<-.,mtxlnG -' 1 ST HF A )JuSTE 
AAA AA = 100.0 
RPP = 1.0 
CCC = 1.0 

ISCALFs * F |. ^ A T < TSI7F-1 ) ) / ( - A XL A T-M TNI. A T ) 

JSCALFs (FLOAT ( JSI 7F-T ) ) / (MfWLNG-MTNLNG) 

51 L = n 

WPIT£<*. 6666)1 A o 
*,4,4.4. FORMAT ( • 1 • .SAfc.//) 

1 L =L ♦ 1 

SIP PfAO<s.)on> xxin, »l »T (i.) , <LON«(L) .jyfap, (X (L.J» »J»1 *4) 
ino FO«-»AT(43 ,fs.?,F«..?.T4.F1.0,F4. i i,Fa.O.FS.O) 

TF ( XL0N«(L ) .1 T.-O0.0. ANO.l, .FO. 1 ) STOP 
TF ( XLONO (L > .LT.-^O. 0) GO TO 40 

WRITE!*. MOO) <XT n **LAT t'L) . »LONG (L) *JYEAR, (X(L.J) « J* 1 . 4 ) 

SJOO F DO * A T ( • • . A 3, ?F7.?. T4.4FT. I ) 

4«| F OP J A T I S X , 4 X , ?F7 . 7 • TS.mFO.4) 

IF IXLATO.) .GT.MAXLAT, om.xlaT (L) .LT.MINLAT.CP.XLONG(L) .GT.MAXLNG 
1 .OO.XLONG(L) .LT.-4INLNG) L=L-1 
TyFaR — JY e 4R 
GO TO 1 
40 CONTINUE 
L*L-1 

WO?TF (6. 13*1 t. 

10*, FORMAT (IX, //.IX, i NU—EP OF STATIONS «FA0 IN = ‘.I10.//1 
L 


0 ^ 324 1 Ml 7 3 1.4 

00 104 I 3 1 , I O 1 7 F 

no 104 j«i • jst/e 

104 P( 1 . J) = 0.0 

on 10 ? *= 1 .LL 

Xl = (HAXLAT-XLAT (M) ) •TSC 4 t. ir *l 

Xj 3 (MAXLNG-XLONG(M) jSCALt'M 

T**I 

IF (XI-FLOAT ( 1 ) .GF. 0 . 5 ) I * I * l 

J 3 XJ 

IF (XJ-FLO.AT! j) .GF.O.S) J*J »1 
lo? P(i.j) s x(-.mi 7 ) • jon.n 

CALL ANAL (IM 7 F. J <: I 7 F,«aaaa.pbp,ccc. 25».01 
00 1021 Ist.TMZE 
00 1021 J 3 1 .JSI 7 E 
10 ?T il( I . J) 3 M ( I , J) - 1 00.0 
J PO l NT 3 ft 
IF ( IPPTNT.FO.) | 

If ALL PONTIN (II. TST 7 F, JST 7 F. 00 .?* 1 .O.LAPI 

C ••••»•••••*••••••••••••••••••••••••••• 

c pf GINNING OF iktfowaTTOn OVER A POLYGON On Tmf map 
PF 4 I NO ? 

PE AO ( ? , l n a ) I^OLY 
S? FONT I NlJF 

c lONLvs THF MIM 4 FP of POLYGONS TO PC PLACED Over The FUNCTION U FOP 
DO ?*« INIIM* 1 , l POLY 
1 nu FOP-MT ( 21 ?) 

s* OF AO ( 2 . 700 | 10 , yl«T (M) .YLONG(M) 

700 FOPMAT ( A* .F*. 2 .FO. ?) 

»I*(MAxlaT-yl«T(M) ) MSCM.f *1 


* P « U 1 


W A 6 0 0 0 1 0 
WAGO 0 O 20 
W AGJOO 30 
WAG 00040 
WAGOOOSO 
WAGO (l Of* 0 
WAG 0 QO 70 
WAG 000 P 0 
WAGOOOPO 
WAG 0 0 1 Oft 
WAGOO l 10 
WAGOO 120 
WAG 00130 
X A G 0 0 1 4 0 
w 46 JO 1 SO 
W AGO 0 1 Oil 
WAG 001 70 
« AGO 0 1 GO 
WAGOMRki 

wagoopoo 

WAG 0021 0 
w AG 00220 
RAG 00230 
WAG 00240 
WAGOOPSO 
WAG 00260 
W 4 G 00 P 7 O 
WAG 002 P 0 
W*GOO?PO 
WAGOO 300 
WAG 0031 0 
GOO 320 


WA 

WAG 00330 
WAG 00340 
WAG 003 S 0 
WAG 00360 
WAG 00370 
WAG 003 O 0 
WAG 003 RO 

W 4 GOO 4 OO 
WAGU04 1 0 
WAG 00420 
W4G00430 
WAG00440 
WAG004rxO 
WAG00460 
WAG004 7 0 
W AG0 04W 0 
4AG004R0 
WAG00S00 
WAGOOS 1 ') 
WAG00S20 
WAGOOS30 
WAG •) 0540 
WAG0OS50 
WAGOOS<>0 
WAG00S70 
WAG00SA0 
W AGO 05 70 
WAG006O0 
waGOOM 0 
WAGO OsPO 
4‘ AGO 0* 30 
W AGO 064 0 
WAG006S0 
WAG00O60 
WAG0ON70 
W AGO 0*P ft 
WAG0OSP0 
WAGOO 700 
WAG007 10 
WAGOO 720 
wAGOO/30 
INTwag 0074ft 
WAG007SO 
WAGOO 7*0 
WAGOO 770 
WAC.O07P.ft 
WAGO07R0 
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* J« (MAXLNrt-YLOWR (M» ) •J^CALF* \ 

!■>! 

J«AJ 

tp m-FLOir u) .fiF.o.m i * i ♦ i 
TF <XJ-FL( iT ( J) .G* .O.s) J*J*1 
KJt.i 

IF (*?.GT. JSI7*-) wRJTF ( 6 • ) 

Ro* F OP** A T (• ? ( >fM T IS OMTSIOF OrtvtFFTIVF F I F.LP IN LONGITUDE •) 

IF (K7.GT. JSI7F) CALL FAIT 

in? su*«u< t ..»> 

IF ( I ,r,T. 1ST 7F) worTF I N.?0<S) 

?g* mo-UT (• ontNT IS OUT^n* 08JECTT Vf FIFLO IN LATITUDE •) 

IF II.(iT.ISI7o» CM.I. f A!T 
TF (SUM.GT.O.) ISMMiSU 4*0.5 
TF (SU^.LT.O. ) JsinsS IM-o.S 
1 0 A T A (HHSIlsP 
IOaTA(Inu b «-»)*Itfao 
IDATM INUM.^L 7 *?) bISU" 
wb l TF (*. ? fc 7) irt.suM 

?A 7 FOOm<\T ( • VALUF AT POINT • « S * , A* t 7 A « • EQUALS » • ?X • F IS. 5 
1) 

CONTINlIF 
3?41 fONTINIIF 

r»0 4444 T«1.TP01Y 

PUNCH 7000 « ( IDATA ( I • J) «J«1 i4) 

7000 F0R-TAT(A4. 14,11,14,14. 15) 

4444 CONTINUF 
GO TO SI 
END 

S«H«0uTTNF «0NTUP(7,M .NJ l.-TN. INT.SC ALK.LA4) 

...PRINTS STATiOAPn MI | MJ (iJT,) WITH CONTOURING OFT4E c U LINES 
...RlN IS THF MjMHtjM VALI'F. INI fs TMf CONTOUOING INTERVAL 
• • • 7£P0 INT mfans ' <n CONTOURS 
...SCAL F IS Tmf SC A(. I N>» FtCTOO FOP PRINTING. 

IF LIJ.GT.?** TwO GRITS AwF PP I jTFO 1 TO ?f> ANO 2* TO NJJ 
OTMFNSINN 17 1 1 00.100) .Z( 1O0. 100* .1 AS<5) 

0 1 MFUS I OM WAI.° ( 16) «L INF ( 1 ?*) *L I‘M?M 

DATA K4LP/1R « 1 p A • 1 H tlHP.lH » 1 HC • 1 H «1 hO«1H • 1 M£ • 1M .1HF.1H , 

1 1MG.1H . I HH/ 

LT0T*INT*1 4 

NJsNJJ 

Jt*l 

IF (NJJ.rtT.7N) NJ«?6 
DO 10 1=1. M 
DO 10 1 • NJJ 

in U(t.J)=7M.J)* c CALE 
IF (INT) SI. 50, SI 
ai NIMsNI-1 
k 0 NJM.NJ- il 
501 PRINT QIO.I.AP 

PI 0 FORMAT ( « 1 • .SA4.//) 

NU5*5*N IV I 

■10? POINT POO. ( 17 (1 .J) • J*J1«NJ) 

O00 FORMAT (?X.?MS) 

00 1 :P«2.‘iT 

oo ? J 0 » 1 .? 

IF ( J1 ,u, n GO TO 70 
00 3 L 8 1 *N J 

3 LIN(L) ■( ( 17 « IP.U-I7(IP-1 *L1 ) • JO) /3* I 7. ( I «- 1 .L ) 

GO TO 30 

?n r>n 40 l = p^«n.ij 

4 0 L-!N(L*?5I>( ( IZ(T0,L)"I7( I 3 *l»L) 1 • JO > /3* I Z ( I R- 1 »L ) 

’0 K»1 

00 4 Jsl.Nju 
LINJ»L IN( J) 

L INF (5) *L I V J 
ND2«LIN(J*1 >-L»NJ 
00 S L * I *4 
K»K • I 

s LTNF(K)=(NOZ*L)/5*LINJ 

Kr5* 1 

4 CONTINUF tn., 

LINF(F).LIN(UJM.1) ^ V;,. 

00 k L«l »N"R Cty* * * -f , 

JOF»L]NF (L)- M !*i P Orir '' 

IF ( JOF ) 4.p,o t'jc 

A JPFbJOF-LT >T* ( (J0F*1 ) /ltot-i ) <4 % '*• \» 

P J« JOF/INT '•n' 

IF(J-1MS,?N.?N r 


WAG00*00 
MAG0051 0 
W4G0O«?0 
WAG00S30 
W AGO 0 *44 0 
4 AGO 0 *4S0 
WAGOOSnJ 
WArtOOM 7 3 
W4GOO«00 
WAGOOARP 

N agoopuu 

nAGOOPI 0 
4AC,0nP?(i 
RAG0OP30 
WAGO0P40 
PAGOOPrO 
WAGOOSNO 
W AGO ON 70 
MDGOOPHO 
WAGOOPPO 
4 AGO 10 0 0 
W AGO 1 C 1 0 
WAG010P0 
WAf.01 P30 
WAGO ) 040 

wagoioso 
WAGO 1 OSQ 
W AGO 1070 
WAGOIOSO 
WAGC10P0 
47G0 l 1 OP 

wagoi i lo 
W4G01120 
WAC-01 1.30 
WAGO) 140 
WAGOI iso 
WAGOI 1 NO 
w AGO 1 ] 70 
WAGOI 150 
WAGOMPO 
WAGO i >00 
WAGO 1710 
WAGO 1^20 
WAGO 1 230 

WAGOI ?40 

WAGO 1 2S0 
WAGO j 2N0 
WAGO 1270 
WAGO 1 ?so 
WAGO 1 ?P0 
w AGO 130 0 
w AGO 1 3 1 J 
WAGOI 3£P 
WAGOI .330 
WAGOI 340 
WAGOI ISO 
WAGOI 3*0 
WAgQ 1370 
WAGOI 3«0 
WAGO 1.390 
WAGO 1 400 
WAGO 1 4 1 0 
WAGO 1420 
w AGO 1 4 30 
4 AGO 14 4 0 
WAGO 1450 
WAGO 14*0 
• AGO 1*70 
WAGO 1 4*0 
WAGOI 4R0 
WAGO ) soo 

wAGolsiO 
w AGO 1 S«0 
WAG01S30 


WAG01SH0 
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S 

*01 

OOl 


7 


404 


S«S 

so 

on? 

007 

S? 

ono 


J*J-(.J/1*> M* 

LISP (H ** AL°( »• n 
o->|NT OOl ♦ ILlNEIL) »L*1 «NUM| 
FO-» *AT I 7 * » 1 2* 4 1 1 
CHIT iMliF 

PPl'iT ono«(T7(l°*Jl*JsJl* Jj) 

CPNTlWlF 

»F (NJ.fO.N.JJ) MfT 'JR* 

Jl«2* 

60 TO SO 
PwpjT 010, LAP 
COUT IMir 
no S? »«i ,\<T 

POINT on?, i 17(1 ,j) ,J»1 »NJ) 

FO->oaT( 'J*.?MS) 
point 901 

TOOT I n ( i ie- 
FOO-UTOh t 

IF (NJ.FO.slJ) octuo- 
NJ*NJJ 


J 1*76 

c,c. T 


TO tfl 

F~n 

SOHOOIITTMF ahai ( 4 I . N ) . A A . AIF 2 • tl_F 4 , MA % ? AS .ERR • Ul. ' I ) 

VARIATIONAL analysis “y «v I T 0 AGnFO pur, nFPlVAllvF FILTERING 
UO IS DATA IJ IS A <A|. YS 1 S TOf.TS 7FCO UO AS *'0 oaTa 
bO'JTINF »It.L AN'AL v 7 r F*)H *OFAS if *.() OAlA 
A6,A|F?,ALF4 43F F1LTC0 wEJ'.mTS BEEE»FN’Ct K.WA6NF4 FDR APPROPRIATE 
VAL'JFS ma.<r AS IS maximum no, of ITERATIONS for is appro* 

C J l TFO I A , o I Nrss I 0 JS OF A-v-AYC ARF IMNl ,NJ) ,U0 (N| ,nj> , ,A (NI ,WJ» 
V a INI •?,N.1«?I WMCAF til 4un N.l APS The INPUT AOWAY S I 7E 
To T S a OUT INF MJMJM17FS T*»F INTEGRAL (T-f SOUAPF.S OF Th£ DIFFERENCES 
I.F Top 


4A601S90 
44601*00 
44601*10 
44601*20 
0Al»O|*lO 
WAC*0l*40 
w AGO 1 oSO 
4 AGO 1 n&O 
WA<j01 S70 
04601**0 
WAoOl NVO 
W AO 0 1 700 
wAGOi 7;o 
04oO}72Q 
44001730 
rf*601 740 
0 AGO l 7S0 
0AoUl 7*0 
0AGO1 770 
0460 i 7PC 
WA6017N0 
4 AGO 1*00 
44601*10 
44G01*?0 


UR 40 I F N T • 4L F ? 

IF T-F L API. AC 1 A“* AI.F4 ) 


«WF *Wf IlyMf < FOP THF ANALYSIS, 

OF -sONTTupF )n, T ypical papraETERS APE 44*100 


Thf Sa'Japfs 
♦ T»-f SOU4PFS 
whFRF aa.alfp.ano alFa 

FO? MFSOSCAIF analysis OF MAGNITUDE 1 
41 F?*l ,0,4I.F4*1 ,0 .RAiUASayp.FPH*. 

N| , ANON J AOF S 1 7F OF APOAY OF UO i O |l 

INCREASING alF > AND 0 3 alf* REDUCES m|GHE« SPEOUENCIES 

Typical max VALUES A&f ALF?* 1 o, o, AN D A|.F** 10,0 

FOR AAalOn.n 

THF VAL"FS OF «l FA? ANO ALE A 4 USUAL 1 Y A ?F ,l«l.»OR 10.0 

DIMENSION U (1 00. 100) «UO( 100, 100) ,YA ( 104« 104) «WA(100« 100) 
FO'J l V Al r MCF (YA ( 1 . 1 ) ,44 ( 1 . 1 ) ) 

10«S 

NjP>*M )•? 

NIP?»NT *? 

Ml*l »N1-1 
*1 1 M ?»N I -7 
N.JM 1 *NJ -1 

NJP 1 *NJ» 1 
NIP l «N 1*1 
PFT A *?,n 

INITIAL I7E GUESS FIELD 4T AVERAGING 

DO IS J« 1 .N.l 

no is l • l . ; ! 
ui i, j» »uni i , j) 
nn io j«i*Njp? 
no in 1 * 1,4 ip? 
ya ( j , J) *o,o 
k Ai r * i 
rONTlNUF 

ChFCS FOP numpfo of no PUPSS 

IF (KNT) IS. 700. is 

KNT *0 

no 12 j* ? . n ip i 
no 12 I«7.»!P1 

YA ( 1 . J) «JI 1-1 . J-1 > 

no op j*?.sjpi 
no op i«?.sTpi 
IFiva u . j) .ftr.i 
sui J *o,n 
C‘.T*n,n 

AVG NINE POINTS 
no P7 JK*1.3 

no 07 is*i,i 

!I«I-?.|K 


is 


in 

?n l 


IS 


1 ? 


os 


24 


.n.nonnn go to ss 


* ago 1 *3o 

4 AGO I 440 

WAGOJP.SO 
4 AGO! AGO 
• AGO 1870 
<• AGO 1 *40 
0 AGO 1 **P0 
4 AGO 1 400 

• 4 AGO 1410 
0AGO 1 920 
4 aGO 1930 
04601 940 

• 0 4 460 1 OSi) 

4AG01 900 
4AG0197Q 
0AGO 1 9*0 
04601990 
4AO02000 
W6G0?O 1 0 
04602020 
44602010 
446020*0 
RAG070SO 
4AG02060 
•4602070 
•AG020MO 
44 G 07090 
0 AGO? 1 00 
0A6O21 10 
0A6O212O 
44G02130 
0 AGO? 1 4 0 
4A6021SO 
046021*0 
4AG0?i 70 
0 AGO? 1 GO 
4AG07100 
0A6O??OO 
HA6O??10 
44002220 
44G07?30 
4AO02?4i i 

4 AG022S0 
4AG0??no 
4AG02270 
446027*0 
4AG0?2 , >0 
4 Ao02300 
44002310 
44602320 
04007310 
44G02340 
4AG023S0 
4AG073NO 
4AG021/0 
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JJ«J-?» JK 

! c CTA ( H , J J) .L c .0.000ftl» GO TO Q7 
os Sll «*>1 )m*va , 1 1 ..)j) 

C*JT«OT« \ . 

07 mntiuhf 

IF ICNT.I F. “>.*) r,o TO o? 

IF (SjM.LF.o.nnoni ) oo to u 

00 TO os 

USF .0001 t^'STrjin OF 7P»0 A VG 

Gi I - 1 * 1 > «.onoi 

oo TO OO 

III 1-1 .J-1 ) *S|.m/CNT 
00 TO oo 


04 

os 

O? 


*s 

oo 

?oo 

I 00 
70 


11 


KNT «<*|T » i 
00 TO oo 

uil-1 . j-i ) < i . j» 

CQrtT I.NiJF 
CONT iN'ic 

wojtE <10. 1 0*') •"iT 

F00-AT(1^0.IS.?«.H points UNSPFCIFIFn This PASS) 

IF |KNT) ?0|,7Q.?0| 

COwT INI IF 

Smooth fitlo of avfpagfs 

00 1| J*?.\JMJ 

00 11 I*?.M-l 

MA(i.j)«< 4 .*Mn.j).tj(|-i, i-ij*u<t*». j-!)*u(!*i.j*])*ijii-i«j»n 
i*?.« mi i . j-n ♦"( I *i. j» ••/( t * j*i i «ti ( \ - 4 . jm/is. 
no i? 

PA 1 1 « Jl «l*-.oii( 1 .Jl *2 .*<i»<?,J) *ij( 1 • J- 1 > *U( J* 1 ) ) *U(2. J-l) *U(?. 
I) ) /IS. 

1? WA (*JI • J) ■ ( a.*i|(NI • J) ♦ ?.• (JlNU'l • J) ♦Ul'Vl . J-l I *UCNI « J* 1 ) ) •MINI*} 
1 ) •OiNTMl , j. ) | ) /IS. 

00 .13 

WA ( I . T) z I . J » ♦ ?.• CJ( I . ?) •<)( 1-1 « 1 > *111 I ♦ 1 , D) *U( 1-1 .21 »Ut 1* 

*»)/1A. 

33 WA I I *NJ) > ( a. T.MJI •?.•!>•( I .NJ'MI «u( 1-1 *NJ> «(M I »1 .Njn »IJ( 1-1 . 
11-UII-l «MJ- ID /is, 

WA ( 1 »NJ) « < 1.*'J( I .NJI ♦?.•(•!(?, mj) *lt| l-NJ-l ) ) ) ) /?. 

WA INI .NJ) « ( 1.*'MNI ,*IJ) • ?.• f • I ( A. J ’1 .No) ♦ * J ( N I • N J- 1 ) ) ♦U(NI >M .NJM1 ) 
WA |N| • 11 s (1.*IJ(N| , I ) •?.• (||(N|M| , I) *UIH! .?1 ) ••MKI1M1 •}) ) /M. 

' * Ml. 21* ~ 


c 

30 

C 

C 


•»S 

so 

I") 

9 


AS 

41 


AS 


WA ( 1 . 1 ) > ( 3 . •* I ( 1 « 1 1 *2.* ( U I ?. 1 ) l C 
00 IS J*1.»J 

00 is 111 , .1 

U I I . J) *WA ( T . j) 

WPITF (T0.3011 

CO NT IM'JF 

FOP *AT ( 1M0.P7H SMOOTH has «Ff»J CALLEO) 

OTMFMSIONAI. analysis 
IN lFP ! OP POINTS AWO wF 10'lTS 
0P0« INI -4 1 * (*' 1-4) 

ALFsR»ALF4*MF TA 

III 10*4. *AI F?.70.*ALF4 

UT 1 J1*-ALF^-m.*6LF4 

ww I T E I TO. Son) AA. ALF?.4LF'..I)I JO.UI 1J1 
ITFOaTi vr sr*.F <f 
0<l 41 ITml.vAFPAS 

1 A« 1 
SilHaO . 0 

00 4? j*3.NJ«o 
00 4? l«3.NT 

U1 J» 01 jo 

C«FC* FOP o°Sf-viTION 

IF I *10 < I . J) .GT.rt.nnooi ) GO TO 43 


) •'J(2.a> ) / m • 


47 

4A 

4? 


AL *0.0 
oil to as 

AL* AA 

uTj«Ul J* Al 

FOUATTON fop pfsioiial 

OFS«-AL*llO I I . J) .1)1 J*0( I . J» ♦nil Jl* lul 1*1 . J) •»!( I-! *JI *ui I* JM ) 
1 - I 1 1 • ALF*Ro III I I - I , J- 1 ) *U 1 T - 1 • J* I ) ♦•! I I • 1 • J- 1 I 4lJ I I • 1 • J* 1 ) I 
0FS*PES.AL« , 4*|H( I .J.?) .Ill 1 ,J-?I «UI I »2.J> *U( 1-2.JI ) 

PL AtPa I ./111 J 

CoPPFC T GIIFSS OF Ij 

1)1 1 • Jl •«!( I . J)-Pl A«P» IPS 

i rMfC* Fno APPPOf j NATION SATISFIED AT ALL POINTS 
F 1 APS(PFS) -FPPI 4S.4S.47 
A*? 

UM*SI|W»PFS*PFS 
CONT INUF 


W AGO? ISO 
WAGO? )«0 
WAG024O0 
WAG0?4]0 
WAGO?4?0 
WAG0?43O 
WAG0?44G 
WAGoJsSO 
*A(.0?460 
WAGO?* 70 
«AG0?4>«n 
WAGO?400 
WAC.O?SOO 

* AI»0?S { J 
W4G0?S?0 

wago?s3o 

waC.0?S40 

WAGO?SS0 

wAG02SsO 

WAG0?S7O 

WAGO?Sfti) 

WAGO?S«0 

WAGO?SOO 

wago?sI o 

WAG0?S?O 
WAG0/S3O 
WAG0/S40 
W4i»0?sso 
J*1 WAGO?SSO 
WAG0?S70 
» J-»WAGO?SHl) 
WAG«T?S«,)0 
WAGO? 7 00 
1.? WiG0?7\0 
WAGO? 7 ?0 
WJM 1 m »i>0?7 30 
WAG02740 

w ai»(»?750 
) /S.WAGO27S0 
WAG02770 
WAG''?7S0 
WAG0?7O0 
WAG02h00 
W AGO?® 1 0 
WAG02M20 
waGO?m 30 
WAG0?*40 
WAG02HS0 
WAGO?MSO 
WAG0PS7O 
W4G0?SAy 

* 600?*lSO 
WAGOPOOO 
WAG02O1 0 
WAG0?Q?0 
WAG02N30 
WAG02040 
WAG02PSJ 
WAG0?OS0 
WAG0?o7 j 
wAGO?OoO 
WAGO?ooo 
WAGO 3000 
waGO 30 1 0 
WAG03O?G 
waGO JO 30 
WAG03040 
w AGO 30S0 

( I • J«4(jO3<i*0 
WAG03070 
WAG030HU 
o AG03PO« 
WAG03100 
w AGO j 1 lu 
WAC.03 1 20 
wAG 031 JO 

WAG0J140 
w AG0 1 1 SO 

* AGO 3 1 SO 
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C WOlTF(p,«,i0| tT.STO 

GO TO (0,* | l«! A 
4 1 COMTIMU* - 
q wPirruo.Ai i) it 
CONT 1*41 tr 

SI 0 F09MAT I 1 *• * 1«F1?.S1 

<11 FO«-*AT I I «, I 7HNO. CF ITFJATIUNS, is» 

S*0 rnw^TlU.M *f |f,HTS/S* K.?/2 *h STO HFVlATION OF rfES!0«)AL ) 
DFTUPN 
F‘J0 


WACOM 70 
WACOM *0 
WACO j 1 s*»» 

wtr,oi/»nn 

V*G03?10 

w*r, 

WAG03230 

WAG03240 

W4C0J>S0 

«AG03?Sft 

WAG0327O 
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: T*IS siirrOUTINF C0*«MTFS TmORaiTMWA I TE COEFFICIENTS 

SURRQUTTnf T*-"3RH (T « I * A . L A T . * 0 J • fJ ET ) 

REAL T ( 1 ?.) , 1 .LA T . A f )J( 1 ?) * J ET ( 12) 

IF ( T ( 7) .GT. 4 0. 0) CALL ICO U T) 

W- 1 1 TF( 4,6*01) 

V* -.TE (6.600?) I AT 

*00) FOrf'iATt • ] • , 'ThoonThwaTTE COEFFICIENTS AS EST I MATEO FROM The 

ll.’FS* > 


12 


I = • «F 1 0 .5 » // ) 
AO J T 


, 0 ) **i 


PET • ♦/) 


GOO? FOR“AT(* 1 « * OF MtAN MONTHLY TF MHE® * TUPF FOh A STATION AT 
1 t.« TnilDF. t .//) 

1 * 0 . 

on s j=i.i? 

TT=T ( J) 

IF(TT.I F.o.jr.n TO S 
Is I* (TT/S, ) ♦ ,S14 

5 CONTINUE 
A1 =0.675 
A2=-77. 1 
A3=l 7Q?o . 

A4*492 IQO. 

As(Al°I**T*A?*»o>?*A3*I»A4) / 10**6 
■WO I IF («,. GOO*! A . T 
GrtOO FORMAT!* ».*A = *fF 10.S.« 

4? I TF (4,6034) 

6004 FORMAT (• •.* HON OL 

on 10 j= 1 • 1 2 
TF(T(J1 .LE.O.»PFT(J)=o. 
lF<T(J),LF.n.)30 TO 6 
PFT(J)=lG.n*(lO.O*T(J)/l*l 
G AOJ( J) sOLHAO <LAT. J) /l?. 

PET(.1>=OFT(.|)»AOJIJ) 
if (T ( J) ,LT. n.)OFT ( J) =0. 

DL=0LPAR (|. AT < J) 

T(.J)=T(J)*l .M.32.0 
PFTf J)=PFT(J>/aS.4 

WRITE (6.6006) J,0L.A0J( J) t T (J) .P£T£ J) 

6006 FOOMATt* • . I3,?X,4F».2) 

10 CONTINUF 
OF TURN 
END 

FUNCTION 0L°AO (LAT. MONTH) 

: COMPUTES average OAYt.ENGTH RY MONTHS FUR GIVEN L AT I TUOE 

OF At. LAT 
INTFGFS .10(11) 

JOUIaO 

JO (?) =3 1 
JO (3) =SQ 
jn ( 4 ) son 
jn cs) =i?o 
JO (6) =151 
JO (7) =101 
JO ( n ) =21? 

JO (9) = ?t 3 
JO ( 1 0 ) s?73 

JOUl ) = 304 
JO ( 1 2 ) =334 
J0(13)=3GS 
JO 1= JO (MONTH) ♦ 1 
JO 2 = JO (month* 1 ) 

OLPARsO. 

on s isjni.jn? 

IF (LAT.RT.4 0. ) DLPA0stlLHAO*0AY4l (L AT. I ) 

IF ( L AT .LF .40 . ) OL R AR*OL a AR*0 A Y4 0 ( L AT » I ) 

6 CONTINUE 

OLflAHsOLRAP/ ( J02-J01 ) 

RETURN 
END 

FUNCTION 04Y41 (W.J) 


C M I00010 
C h I 0 0 0 2 0 
Cm I 00030 
CHI00040 
C* 1 1 00 050 
C*1 00060 
V ALCN 100070 
CmI 00020 


O'* 

0/^ir 

°0#^ P 4 






• »F5.2» 'NCMIOOORu 
CM I 00 100 
CMJ001 10 
C M I 00 1 20 

cm i oo no 

CM100140 
Cm 1 00 1 50 
CMT OOloO 
C M I 0 0 1 70 
CMIOOIMO 
Cm 1 00 1 RQ 
CMI00200 
CMI00210 
CMI 00220 
CMI0O230 
CMI 00240 
C M 1 00250 
CMI002GO 
CMI00270 

cmiooppo 

CMI002SO 
CM 1 00 300 
CM 1 00 310 
CMI00320 
CM I 00 330 
CM I 00 340 
CM I 00 350 
CM 100 36o 
Cm I 00 370 
CMI003OO 
Cm 1 00 3OQ 
CMI00400 
CMI00410 
CM I 00420 
C M I 00430 
CMI 00440 
CM 1 00450 
CMI004GO 
CMI 00470 
CMI004P0 
CMI00400 
CMI 00500 
CMI00S10 
CMI 00520 
CM I 00530 
CMI00S40 
CHI 00550 
Cm I 00560 
CmI 00570 
C M I 005 p d 
CMT005Q0 
v.Ml 00600 

cmioogio 

CMI 00620 
Cm I 006 30 
Cm I u 0 g 4 * j 
C*M00650 
CM I 00660 
CMI 00670 
CMI006H0 
Cm I 00690 
CmI 00700 
CMI00710 


TmIs FUNCTION SuRUOUTInF COhRMTFS OAYLENGTHS FOR 
L ATT TUOF S NOPTH OF uf) OF GHFF S • 

04 Y4 1 *12. 25* (1 .6164* l . 764 3* (TAN (3. 1 4 16*»/180. ) ) * # 2 > *C0S ( 0 . 0 1 72* J-2 Cm I 0 0 720 
1.R4) C M 100730 

RETURN CMI 00740 

FND CM I 00 760 

FUNCTION I)AY40(W.J) Cm I 00 760 

CMI00 »’ 70 

THIS FUNCTION sumoouTINF COMPUTES OAYLENGTHS FOR CM 1 0 0 7 8 0 

LATITlJOFS AT OR SOUTH of 40 DEGREES. CMI00790 
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FILM CmISURS FORTRAN A PUPniJE / L ARS .1031 


A. 
«. 
r. 
o. 

F. 
F. 

£ o. 

C m. 
C 1. 

C vJ. 

$ f: 

r **. 

R. 

C P. 


10 

IS 

>0 

?s 


10 


IS 

AO 

AS 


n»YAO* 1 ?. 1 A* I 1. 1!«*T 4N» J. lA)S«*/ 1A0. ) ) *COS <0.01 72* J-2.R4) 

or turn 
F\n 

FUNCTtiV pp T (T.T | ,Ta, T f 1 

C0 -o lJT c S T-ORT TmkaITE 'S ESTIMATE OF POTENT I AL F V APOTRA nSO H AT ION 
IF ( T «LF . <i • ) PF T:0. 

IF < T .I.F .0. ) he TURN 

PFT = iH.o*rr*< no,n*Ti/Tn»*TA 

return 

EnO 

SIJRIOUTTA'F TfOM(T) 

PF AL TM21 
no s is i.i? 

T < I ) s (T <1 ) -1?.0) /1 .H 

pr TURN 

FKP 

SlJRROUT I NF Hvnw« < T.o.SPS.SOU.SP.^.PP.PF.PL.FT .W.L«Hl . TA.TC) 

DmS PaLmEpiF myOOOI OSIC *.« COUNT IMP f i'iw ONE MONTH. 

T: The mean MONTHLY T-*Hpr6« 1 U->f T ►. OF t.PFFS F. 

p: The total monthly ppF r i p i ta r ion in inches. 

SPS: initial Amount of ayailahlf MOISTUPF 

The SOU IN INChfs. 

I N l T I A » amount OF AVATLfRIF POISTUPF in the 
OF T«F coil IN f NC mF s . 

INITIAL AMOijnT OF AvilLAPLt MOISTURE IN ROTH 
SOR IN INCHES 0« <sns*SO|l. 

The avaIlaplf moisture capacity in inChfc minus 
potent !A(. oe rnsOi-E In lNC Mt 'S. 

POTENTIAL FVAP0TOAMS°I u AT ION in INCHES. 

POTENTUl I OSS OE SOIL I»r\ I ** T • I s-r F | ri IA'ChES. 

COPPUTFO Fw APOTP ANSPI JiTIvt , in Inches. 

PprM aU.-.F 5 net SAIN IN SOI! MOISTijRF IN INCHES. 

NET I OSS nr SOIL moISTu-'F Jn InChFS. 

T mopm Th* SITE’S hF i T [nof*. 

THOPVTHMAI TF is COEFFICIENT COMPUTER FPOM Hi. 

ThOPNThnaI TF »S MONTHLY AOJUSTmENT eon 0AYLEN3TH. 

T «b»SJ,PW,ue 


SPu: 

S°: 


.B.l .LS.LU 


Ta.TE) /?S.a 


1) 


w: 

Drj ; 

pe : 

pi : 

FT : 

R : 

L : 
hi : 

Ta: 

TC: 

real 

peal PL. ft 
S paSPS»SPU 
PRsM* 1 . 0-SB 
TT * ( T-l?, 0 ) / 1 
PF=PFT < TT.hI 
A=SPS-PF 

IF < A)S. 10. 10 
P| SsSPR 

PLUS (PF-PLS) *SUll/ ( *• 1 . 
if <SP||,LT.PLU»PLU = SPU 
so TO 1 S 
PLSsPE 
PLU*0 . 

PL=PLS*PLU 

P=p-PE 

I p ( a I ? A , * C • A J 

40 r A MS <0) 

IF(SPS-AM) js.30,30 

LSsSPS 

BSsn. 

SS = 0„ 

LUs ( AO-LS) •SPo/ (M* 1.0) 
R'lrO. 

Si'sSPij-LU 
SO TO 3S 
L S= AO 
PSsO. 

SSsSPS-LS 

LU*0. 

RllsO . 

SUmSOU 

POr.1, 

FT=h*LS*LU 
GO TO ?s 

C*l .0-SPS 
IF <C-M) *S.SS.*e, 

BC*C 

LS«0. 

SSsl .0 
OsH-RS 
Fstf-SPU 


IN THE SURFACE LAYFA 
UNOFRLV IN6 
LAYERS OF ThE 
1.0 INCH. 


C M I OOMOO 
CMIOIUUO 
CMI 00A?0 
C m 1 0 0 R 30 
C ’* I 0 0 h a o 
Cm I U0R80 
Cmi oohnO 
CmI 00870 
CMI OOH«0 
CHI 008R0 
Cmi OOMOn 
CMI 00R 10 
CMI 00 l?0 
CP I OOP 30 
C*l nnR*o 
CMI oooso 
C^I oorho 
CM l O0R70 
CMI OORMO 
CHI00RR0 
CmI 01000 
OFCMI 0 1 0 1 0 
CMIOIOPO 
LaTE°Cm I 0 1 030 
C 1 * 1 01040 
C“I 010SO 
CMI 010O0 
CMioioTO 

v'Ml 01080 
C"I 01 OM0 
CNI 0 1 ) 00 

CMfOU 10 
CM 1 0 1 IPO 
CM I 01 1 30 
Cm I -J 11*0 
Cm I n 1 1 SO 
Cm I o 1 1 on 
CMI l) l 1 70 
rMioi inn 
CMlOllRO 

c*i oi?oo 
CMI 01210 
Cm I 0 1 220 
C M I 0 1 230 
C«1 012*0 
Cm I 0 1 2S0 
CMI 01280 
CM I 0 1 2 70 
CMI 01280 
C“I 012 Q 0 
CM 01300 
Cm 01310 
CMI 01 320 
Cm 01330 
CMloi 3A0 
C- 013S0 
Cm 01380 
CM 01370 
CM 013M0 
C M OIJRO 
CMIOUOO 
CM10U10 
CMI 0 1 a ?0 
C M 101*30 
CMI01440 
CMlOUSO 
Cm I 0 1 *60 
CMI 01a 70 
CmI01*m0 
CMI OlARO 
CMJOISOO 

cmioisio 

CM01820 
C m | II IS30 
Cm l 0 1 S* 0 
C M 1 0 1 SSO 
CMJ OlSnO 
C M 1 0 1 S 7 0 
CM101SR0 
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IF (E-O)SO.SS.SS 
cr) omiF 

U0*6-R'l 
r.o TO *g 
«S WllaO 

on*'*. 

*o t.u»n* 

Sl'«S*»U«PU 
GO TO 70 
AS PSs* 

LSaO. 

SS*^P5*°S 

Mils'). 

L*i«0. 

SiJaSPU 

POaO. 

70 FTaPF 
7S PrRN*Rn 
l*LS*L'i 
SOSaSS 
SPU=SU 
SP=SPS*SPU 
RFTHHN 
FMD 

SIJHROlIT INF MAD” (Y . YP 4R. A.IR.PhaSF ) 

INTEGE 0 T.) / 1 ?/ 

°FAl Y (1?) . AMP (A) .PHASK (6) «»(M 
DO >1) Kal.6 
Si» M A = n. 

SUMMcn. 

«;iimyso, 

SuxYSu»n . 
p TO. 1*1* 
no S Tr I ,L 

«;iiMA2Si!'<A*v(T)*rr>S('<o?»Pi<»T/L) 

Sl»MRaSUM Q *7 (T) »<; INI I* T/L) 

SU'*Y a SU *Y ♦ Y ( T ) 

Sil”YSO = SU»- rSO*Y(T)**2 
s continue 

A = 2*SUMA/I. 

0 = ?*SU*'D/i. 

YM AW = Siimy/L 

YVAOa (SHMYSP-I • YM AM** ? ) / ( |_ - l ) 

amp (< ) aSOOT ( a»*?.p»*2) 

Al PmA=AT*N I APS t A/P) » 

Phase <k ) =al jm * 

IF ( A.OF.O. ANO.P.LT.')) PHASP (K) aPI-ALPHA 

1 F ( A .L T.o. ANO.-..1 T .0) MM ASP (K) sPl ♦ ALPHA 
IF ( A.L T . 0. ANO.u.OF ..)) DMA Sr ( K ) a?<*P I - ALPM A 
P|K) a < A*P(< )••<?)/(?. 0*YVAP) 

ARMASFrPHASF <" ) • 1 “O./PI 

IF ( APHASF. 3F. ) on.) ADHASFcAPMACE-3fO. 

AMGLPs^n./l 
SOYaSOPT (Y'/AD) 

WUITF (A.M)0rt) < 

6000 FORMAT t • •. ‘STATISTICS OF MaPMOnIC NlIMflFO: • • I ?»//) 
WPITF (A.600S) 

aoos format (• ‘.‘mfan OF y amP(K) angle phase 
1/) 

WRITE (A. 60 10) YRAP.AMP(S) . ANGIE.APMaSF.P(K) .SOY 
6019 FOP- a T I • • .2F9.A.FA. 2. F0.2.F7.*. FI 0. *•//✓//> 

70 CONTlNMF 
PFTUPN 
ENO 

SIIPPOUTTNF CM | (LAT.W.mI .T>. ARAM. A A”P. APHASE) 

C A. COmPUTFR *FF*(_Y VAL'IFS OF C“l. 

C °. iaT: iaTITuof in ofgdffs imunowftmS) . 

c C. w: avatlam L c watfr CAPCITY in INOES minus 1.0 INCH. 

CO. hi? tmORntm* a I TE »s Hf aT inhe*. 
r f, m: thoontmwa i tff cofffic.ifni computed fro** hi. 

c F. AO A 4 S THF MEfi.| VALUE OF AI.PHA FoD ALL “ONTmS. 
c G. A AMO : THF A”t»LTTUOF of ThF K Tm HARMONIC OF SL°ha. 

C h. aphasF: the phase angle OF Thf * th hapmomc OF alpha. 

nATA PT/3.1*|Su/ 

PFAL L»T.I\AMP(A1 .APmASFI*) ,M 

InTfgeo luc is) 

A W (* = W ♦ 1 . 0 

PE A0( A.4 0nn.FNOa'>S)LOC 
*onn format <sa*) 


SO OF Y* . 


Cm I n 1 SQO 
C”IOIM)0 
CMI 01610 
CMI 01620 
CM1 01630 
CHlOlNAO 


CM101NSO 
0*1 0 lt»60 
0-101670 
chi oisrto 
CMI OlhPO 
CMI 01 TOO 
CMI 01 710 
CM1017P0 
Cm 01730 
Cm 017*0 
CM 0 1 7 S 0 
0*1017*0 
CM 01770 
CM 017SO 
CM 0 1 7 VO 
CM 01800 
CM101»10 
Cm 0 1*20 
CM 01*30 


CmI 018*0 
CMI OISSo 


CMJ 01060 
CM I 0 1 S 70 
Cm 101 RhO 
CMI 01 «P0 
C M 1 01R0O 
CM I 01? li> 
CMlOls^O 
C M I 0 1 P30 

cmiojpao 

C M I 01VS0 
0*1 01P60 
CMI01P70 
CMI 01RK0 
CM I 0 1 RQO 
CMI 02000 

cmio?oio 

CMI o^o?0 
CMI 02030 
CMI 020*0 
C M I 02-3SO 
Cm I 02060 
CMI 02070 
C M I 02080 
Cm 1 02 OP 0 
CM I 02 1 00 
CMI 021 10 
CMI 02120 
CMI021 30 
CMI 021*0 
C* 1 02 1 Su 
0*102160 
CMI021 70 

CM } 02 l SO 
O* I 02 1 VO 
0*102200 
CM I 0?? 1 0 
0*1 02220 
0*102230 
Cm I 022*0 
CM1022S0 
C“I02?6i) 
CM! 02270 
Cm! 02 ? 80 
CMIO??<JO 
CMI 02300 
CM1O2310 
C * 1 02 320 
C‘M 02330 
CMl 023*0 
0*1 02 3So 
Cm I 02 360 

C M 1 w? 3 ?o 
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C A. 

C ». 

c r. 
c n. 

♦ no i 

♦♦no 

*001 


son* 

s 

♦ nos 


,//) 


in 


PF4DU.4 00) * r fc li =OF)SPS.SP i.yp.gm 

SPS: in j T I al sh-FaCE son moTSTiiof CONDITIONS in inches 

SP'J: iNtllAl SUBSURFACE SO I l t «iT'l»*F CONDITIONS IN INCHES 

vP: initial VALUE OF funFx OP t vapOTwansp |paT IOn OEFICIT. 

pp; INITIAL value OF Inop* OF tXCFSSlVF MOISTURE. 

FORMAT (F4.2. acs.^i 

vtP T TF ( *..*nooji or .laT 

Format l • l • • ' wFF"LY Cm' VAL'IFS pop: «,SA4.» LATITUDE: 'tFS.^* 

wo i tf (g.goo 1 ) a.t ,s»s.sPu« yp.gp 

FOPMAT(» •.•INITIAL CONDITIONS: »#//• • •••ARC: ••FS.2»/»' *t*SPS 

1 •.fs,?./.» • • • SP'J : ».ff.2»/.» • . • vm: • •F‘>.2 # /.» '.‘GP: 

1P.//J 

WOITF (6, 60 OF) 

FOP-ATC •.' wF T P OF PE ET ALPHA CET VP rf SP 
1 SP C M I'./I 

PF ( 4 . 4 0 O^ .F'H'sQS ) JO. T .P 

FOP '•AT H3.F3.n,4X.F4.1) 

IF (L AT.0T.40) TF=7.0<»DAV4l (LAT.JO) /3* S.2S 
IF(LAT.LF.4O)TC= T .0»OftY4')(LAT,JP)/-»SS.^S 
CALL HVOPA ( T .P.SOF .S»J.SP.w.PD.PF,PL .FT .H.L« hI.TA.TC) 

A=1?.0<> ( JO* 1S.0) /36S.2S 

on io t=1 . *• 

ALPHAS ALPHA* A A”p ( I ) •SIN(?*PI<*A/ia.0*4PMASF ( I ) ) 

|f ( alpha. or. i ,r>) ai pha 


C M I 023*0 
C M I 0?3O0 
C • ! 02-00 
c**l 0 rf <• 1 n 

oio?*?') 

C*I024 TO 

C M i 02 --n 

CMI O^ASO 

0*1024*0 

C*I02470 
• «F5.C'«I 024* n 
C*f 024*0 
CM02500 
C*I0?S10 
CHioaspo 
CMIOPSJO 
CHI02S40 
Cn I 02S50 
Cm I 0 2S*0 
CMI 02S70 
CMI0/S*0 
C*I 02GR0 
CMI 02600 


)= 1 .« 

If ( ALDHA.lt .0.0) ALPHA =0,0 
OF T = AI Pha*hF 
0F= (FT-CFT) /SOP T r AL J HA) 

YP=n.f.7«yp*i.A<»nF 
M: (S°*SPS*S°") / { ?* <SP*PR) ) 

IF (YP.LT.o. ) v=yo 
IF(yP,GF.O.)Y=M*YP 
IF (GP.FO.n. ) HaO. 

IF(r,D.r,T.O.)Hir.w 

IF (GP.0T.0.5)m=0.S 
IF(r,P.f,T. 1 . 0 )H = n.cor,p 
G=GP-H*m®R*PO 
C = Y*G 
G e =S 
I T = T 

^PITE (S.S010) JO. I T.o.OE»PP .FT . ALPHA «CET. YP.H.SPf GPfC 
*,010 FORMAT (• '. 2TT. FS.l.F*. 2. F*. 2. Fa. 2, F6. 3 «FS. 2, F*.2« 3FS.2.F5. 1 ) 

GO TO =« 

"AS wPI TE ( *,. 60 i S ) 

A01S FOPMAT ( • •.///.• '.'wF: JULIAN DATE OF REEK ENDING. •./«• «.»T 

1 MCJM TFMpFOft T'JPE FOP Thf *FFv (n*G. F. )•*/«» ».»P: TOTACVI 02*20 

1L ppfc t PI T a T ION p op Tmf «*cf< in inches. ••/•• '••OE: RELATIVE TCMin?M3<'> 

IRanfpipaTIon anomaly fop thf wFF*«./.« •.•pf: POTENTIAL Ev apoTC*! 02*ao 

1RANSPIP1TI0-I F oo T-*F .FFK .•,/,• • . • F T : CO'^PHTEO FvapOTpA'IShIPaC * l 02**0 

1 T l o.'i FOP ThF wFF k .•«/.• •.•alpha: FSTIHATEO COEFFICIENT OF EVAFOTCmI q:>m*0 


CMI,>?61 n 
C*! 02*20 
CMI02630 
CMIOP640 
C*I 026S0 
Cm I n 266 0 
CMI02670 
CM I 02*Ai) 
C M I 026R0 
CMI 02700 
CM I 027 1 o 
CMI 02720 
CMI 02730 
CmI 0?7uf) 
CMI 02750 
CMI 02760 
CMI 02770 
C*I 027B0 

Cnl 02700 
Cm I 02*0 0 
Cn 1 02*10 


IQf \ 3 P 1 OAT TOM FOP ThF vprK.'./.' • • • CFT 

wp I TF (*.*0^0) 

60PO FORMAT (• •.•mo; INOFX 

l: average percent of 

ISP: IIMfMfS OF WATER I 

1SSIVF soil MOTSTURE. • */• • 

] THF PFE K • • ) 

RETURN 
FNO 


ALPHA X PF • ) 


CmI 02470 
CMI 02**0 

OF FV APOTPfiNSPIPAT IO») OFF I C I T . • . / • * • . * *C * I 02M9O 
field capacity DUPING The •EE<. , »/«' •. , CmI02m 00 
i ThF foil.*./.* *.*GP: INDEX OF £*CEC M 1 02*10 

'.•CMI: CROP MOISTURE INDEX VALUE FOPCmI02h2) 

C m I 0 2 v J 0 


C**I 02040 
CMI02PSO 


JRIGINAL page ifc 
OF POOR QUALITY 
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FUNCTION OAYl.R(V.J) BOROOOIO 

COMPUTE nfYLFNMn Fop a GIvFN JULIAN 04 TE AND L A T I TUOE 90R00020 

f)AYLG*l?.2 e * (1 .*16 m* 1 . 764 *»(TAN(3. 1 4 1 6*w/ 1 An . » ) • ) *C0S ( 0 . n 1 ;?• J-2 k('RP00 30 

1.04) HOH00O40 

IF<w.GT.40.n)PETURN PO800O50 

OAYLG« I ?• 1 4* ( 1.37»TAN(3. l4l(S*V/lflO. ) ) *COS(O.017?»J-2.94) MOHOO OH* 

WFTURN ROHllOO/O 

FND POROOObO 

SURPOUT1NF S^BMTStPn.LAT.wXt PQH00090 

COMPUTFS Haily SSPWIN6 WMFAT STAGE oaTES FROM OARY TEMREPATUPFSROROU 1 00 


PFAL 


ROMOOl 10 
90HOO 120 
MOHOO 1 30 
90*30 0 1 *0 
MOR0O150 
POHOOl t>0 
MOH001 70 
MOHOO 1 A 0 
ROROO 1 90 
ROH00200 
MOHOO? 1 0 
HOH00220 


POOO 


10 


IS 


?S 


w< (*. A46> 

A*-'° AY MX FONT A T f . ■> ONF YFA9 OF OARY mFATMFR VAPlARLES 
MX ( 1 .N) = OA It Y MAXIMAL: T E k ‘H fc ‘ o a t UPE IN PEG. F. 

WX ( ?• N ) *0 A R Y TFMHfcRATilPF IN RFC,. F. 

MX (4 .Nt tCOVBMTFO BPp I imp WHFAT STAGE 
TNTFGFR PH 

PFAL LAT.4An(«;),PAi(B»,MAP(B) .PRO (9) »RP 1 (5) ,PR2(5) .PCI (5) .PC2(5) 

PFAI. FR TTA(F) • C u I T»* ( 5 ) • CP ! TC ( S ) 

Oil 10 * = 1.S 
CPI TA ( I I sQQOOB.Q 
CPI TH ( It =09030.0 
CPI TC ( It =oqoqo.<5 

PFAn(??,??nO0 1 bao ( I ) .RA1 ( 1 1 .PA2 ( I ) .PP(i ( I ) ,RR1 ( 1 1 .9B2< I ) .PCI 1 1 > .PC2RGH00 230 
l(it PU800240 

FORMAT ( RE 1 0 .4 | HOROOPSO 

IF (OA2( It . wE.O. tCOITA ( T)»-(PA] ( I )-?*RA?( I ) *PA0 ( It > / (?*PA?(T t ) MCM00260 

IF (PR?( T t . wF.p. t COITF. C 1 1 =-( JM1 ( T I t •pan ( I > t / ( ?<*9R2< 1 I > ROR00P70 

IF (PC2(I t . mF. o.l CR I TC ( f ) =- (PCI ( I »-?*MC2( I )«PriO (It) ✓ (?*WC2( I > ) POP 0028 0 

PFAOS Thf FWF SFTS if MO^FMTSOM COEFFICIENTS FOR each STAGE POH00890 
ANn COmouTFS CRITICAL VALUtS OF FACh QUADRATIC POROO'XOO 

CONTINUE POB00310 

F=1 MOR00320 

CtlMsO.O ROHO0330 

010=0.0 POB00340 

OO ?S T=1.36S POPOU3SO 

0L=0AYLG (L AT. T 1 ROROP36U 

IF ( ! .FQ.pnt riiM= ) .0 ROH00370 

IF (1 .LT. (Pi* 1 ) t GO Tn 2S 90»3003b0 

if (r'jM.sF.4.oi n in=o.o Riteoonoo 

I F t CUR • GF • • 0 ) GO TO 2 C R0«00400 

V 1 sOUAO (OL .PA0(K1 .RAt <) • 4 A2 (f ) . C° I T A ( K t t ROROOalO 

IF (PAOtKt .FO.ROQ. J Yl “ . .0 ROR00420 

V2=QiJAn (MX ( 1. 1 1 ,o°0(K> .MSI (Kt ,m 3 2 < K > »C B I TP ( K ) ) HUP 004 30 

VXsOtJAO (W-x (2. 1 1 .PRO (<) .RC1 (* ) .RC2 ( H. ) .CPI TC <K) > POM 00440 

OIO*V1* (V?*wl) RORO0450 

ACUM=C'lM.niO ROH00460 

IF ( ACU (.GT.ft.OtGO TO IS POH00470 

IF ( ACUH.LT. (f ♦ t t t 00 TO 15 POH004S0 

IF(Oln.LF.O.OtGn TO 25 ROR00490 

OIFr (K. t ) * 1 ,0-Cl‘ M POROOSOO 

PCT=0IF/0I3 HOBOOSn 

RRCT= 1 . P— PfT MOH00520 

V] =0IJAP (OL .PAO (K*lt.»Al(*<.lt,PA2(K.l),CRITA(K*lJ» POP 005 30 

IF (PAn (**1 ) .EQ. 090. IV 1 = 1. A MOR 0 0540 

V2=0U AO ( w X ( 1 . T) ,NC a (* ♦ l > ,WR 1 (K ♦ 1 t .RP2 (*♦ 1 ) «CH T T« (K ♦ 1 1 ) POP 0 OSS 0 

V3*CUftn(wx (2. T » .M°0 (K*l) .RC1 (K* 1 » .«C2 (K* l ) tCRITC (K*l ) > W0«005ti0 

0T0N.V1* (V?»VA) ROH00570 

nin=Rci # nn*ppcT»oiftN poboosho 

CUR=CtlR»0(0 * POR00S90 

I F ( CUM .OF, (K.J t )F=K»1 POPOO^OO 

IF(CU“.GF.4.0tCUM=4.l) ROHOOMO 

wX(4.It=CUM POPOO'520 

PFTUPN PdHOOO.'lO 

FNO POR00h<*0 

FUNCTION QUAO ( X . AA, A 1 . A?,r ) 90000*50 

COMPUTFS VAIUF OF PjPfwrSnN'S OU AflP AT I C TERMS. POP 00660 

01 1 A0= A 1 • ( X-AO J ♦AP* (X-Alt »•? HUH00670 

IF ( A2.GT.0. ANO. X.LT.C) QUAOsO. ROHOO6A0 

IF (OUAO.LT .0. t OUAn=0. 9 OHO 0690 

PFTURN MOP00700 

ENO ROR00710 
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APPENDIX B 


The structure for 
algorithms used in 


INTERACTIONS OF SUPPORTING SOFTWARE 

interactions of the supporting software for the models and 
these procedures is depicted in figure B-l. 


B-l 


INPUTS 



B-2 




Figure B-l.- Functional flowchart for the supporting software. 



